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ABSTRACT 



Rotation has a number of important effects on the evolution of stars. Apart from 
structural changes because of the centrifugal force, turbulent mixing and meridional 
circulation caused by rotation can dramatically affect a star's chemical evolution. 
This leads to changes in the surface temperature and luminosity as well as modifying 
its lifetime. Observationally rotation decreases the surface gravity, causes enhanced 
mass loss and leads to surface abundance anomalies of various chemical isotopes. The 
replication of these physical effects with simple stellar evolution models is very difficult 
and has resulted in the use of numerous different formulations to describe the physics. 
Using stellar evolution calculations based on several physical models we discuss the 
features of the resulting simulated stellar populations which can help to distinguish 
between the models. 

Key words: stars:evolution, stars:general, stars:rotation, stars: abundances, 
stars: chemically peculiar 



1 INTRODUCTION 

The effect of rotation on the internal p hysics of stars has 
been considered for many years (e.g. iKippenhahn et al] 
Il970l ). Rotation causes sig nificant changes in t he hydro- 
static balance of the star (|Endal fc Sofial Il978l ), t hermal 
imba lance causes a meridional circulation current (|Sweetl 
ll95Cfj and differential rota tion leads to shear instabilities 
(e.g. ISpiegel fc Zahiil 119701 ) . These all result in the mixing 
of angular momentum and chemical elements within the 
star leading to changes its surface properties such as the 
surface gravity, temperature, luminosity and chemical com- 
position. Over the course of several decades, the physical 
formulations used to describe stellar rotation have prolifcr- 
ated (Zahn 1992; Talon et al. 1997; Mcynct & Macdcr 1997; 
iHeger et al.ll200d : IMaeder fc Mevnetll2 005). Whilst each new 
model has been suitably justified physically, there has been 
little observational data to back up claims of improved phys- 
ical agreement. This leads to the possibility that any num- 
ber of physical models can be chosen to produce a range 
of desired results which may or may not be accurate. This 
situation is worsened because the data required to constrain 
the models is scarce. However, with the observations of the 
VLT- FLAMES survey of massive stars (lEvans et al." '2005', 
I2OO6I) and VLT-FLAMES Tarantula survey (,Evans..2011. ) it 
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is now becoming possible to make such comparisons of dif- 
ferent physical models and place some constraints on the 
formulations used. 

Comparing stellar models is still problematic because of 
the difficulty of isolating the effects of rotation from other 
physical and numeri c al diff erences in the results of other 
groups. Potte r et all l|2012l . hereinafter referred to as Pa- 
per 1) presented ROSE, a code capable of performing stellar 
evolution calculations with a number of different models of 
stellar rotation, eliminating any differences owing to other 
numerical or physical effects between different codes. 

In Paper 1 we considered the main differences be- 
tween the evolution of individual stars under the assump- 
tions of several popular models. In this paper we combine 
that analysis with t he stellar population code, STARMAKER 
(IBrott et al. 

I l2011bl ). to determine the difference in stel- 
lar population s that arise f rom t wo physical models. One 
is based upon iHeger et al.l (|2000l ) and has solely diffusive 
transport o f angu lar m omentum. The other is based on 
iTalon et all l| 19971 ) and iMaedeij l|2003l ) and has both dif- 
fusive and advective transport of angular momentum. The 
two different models have very different diffusion coefficients 
and there are marked differences in the results for individ- 
ual stars. It is possible to get better agreement between the 
models under different criteria by adjusting the associated 
unknown constants but this leads to poorer agreement else- 
where. 
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One particular consequence of difTerent input physics 
that we found in Paper 1 is that the mass dependence of the 
mixing is very different in each case. The models agree in 
that the total enrichment in low-mass stars (M < 20 M0) is 
much less than in high- mass stars (M > 20 M©). However, 
the enrichment found with each model is very different for 
low-mass stars despite reasonable agreement for high-mass 
stars. In Paper 1 we also concluded that the difference be- 
tween the two models varies for different met allicitics. For 
Z = 0.001, the model based on iHeger et all (i2000) actu- 
ally produces significantly more nitrogen enrichment in high- 
mass stars, particularly for slow and moderate rotators. We 
shall explore all of these features further in this paper. 

In section 2 we briefly outline the details of the two 
codes, ROSE and starmak e r. Fo r full descriptions we refe r 
the reader to iPotter et"ai] (|2012l ') and lBrott et all (|2011bf ). 
We also describe the models under comparison. In section 3 
we compare the stellar population predictions and consider 
the similarities and differences between the models. In sec- 
tion 4 we present our summary and conclusions. 



2 INPUT PHYSICS 

Let us present the physical ingredients of the numerical evo- 
lution code, ROSE, and the population synthesis code, star- 
maker, as well as the two models under comparison. 



2.1 RoSE 

ROSE is based on the Cambridge stellar evolution code, 
STARS , the first version of which was written by Egslcton 
(|l97ll ). It has been modified and had its physics updated 
many times since. For deta ils of the last major update see 
IStancliffe fc Eldridgd (|2009h . The code solves the four struc- 
ture equations, seven chemical equations and now the an- 
gular velocity equation in single, implicit, Newton-Raphson 
iterative steps. We calculate ^H, ^He, ^He, ^^C, "N, ^^O 
and ^°Ne implicitly and 39 other isotopic abundances can 
be calculated explicitly. The structure equations including 
rotation and the implementation of the angular momentum 
evolution are described in Paper 1. We summarize the key 
components here. The equations for hydrostatic an d ther - 
mal equilibrium arc modified as by lEndal Sofial (|l978l ') 
and .Mevnet fc Macdcr (199'i^. A surface of constant pres- 
sure, P is defined as Sp, while Vp is the volume contained 
within Sp and rp is the radius of a sphere with volume 



Vp = 47rrp/3. Mass conservation gi' 



dmp ^ 2 

— = A-Krpp, 

drp 



(1) 



where mp is the mass enclosed within Sp and p is the den- 
sity on the isobar which is assumed to be uniform on Sp 
even when the star rotates differentially. Strong horizontal 
turbulence keeps relevant physical quantities uniform along 
isobars. The local gravity vector is 



ffeff 



Gmp 2 .2 
5 V^l rp sm 



Er + rp sin 9 cos 6\ et 



(2) 



where fl is the local angular velocity. The average of a quan- 
tity q over Sp is defined as 



<q>= j qda, 



(3) 



where da is a surface element of Sp. The equation for hy- 
drostatic equilibrium is 



dP 
dmp 



Gmp 



where 



fp = 



4nrp 



< Soft > 



GmpSp 

and gcff = IffcffI- The thermal equilibrium equation is 



dlnT 



3kPLp fr 



dlnP IGTvacGmpT* fp' 



(4) 



(5) 



(6) 



where Lp is the total energy flux through Sp, P is the pres- 
sure, T is the temperature, k is the opacity, a is the radiation 
constant, c is the speed of light, G is the gravitational con- 
stant and 



f A-rrrj, 



(< ffoft >< 



gcfi 



>y 



(7) 



Again, the non-rotating equation for stellar evolution has 
been preserved except for the multiplication by fr/fp- Of 
the two factors, fp deviates further from unity for a given ro- 
tation than fr- Additional secondary effects of the reduced 
gravity must be taken into account when calculating quan- 
tities such as the pressure scale height and Brunt- Vaisala 
frequency. For the remainder of this paper we drop the sub- 
script P on the variables. 

We use the mass-loss rates o f IVink et al.l l|200lf ) for mas- 
sive stars although several other mass-loss rates are often 
preferred. The enhanced mass-loss rate resulting from rota- 
tion is given by 



M = Afn=o 



GM 



r2 < geff > 



€/2 



(8) 



wh ere we take ^ = 0.45. This is the same as the one used 
bv iLangeil l| 19981 ) except that we have chosen the critical 
rotation rate such that QctH corresponds to < g^ff >= 0. 
For a more complete descript ion of the critical rotation rate 
of stars we re f er the reader to lMaeder fc Mevnetj (|2000l ) and 
iGeorgv et all (|201ll '). 

The thermal flux, F oc gcfi{9) (|Von Zeipell Il924l ) 
strongly depends on co-latitude. This produces a thermal 
imbalance that drives a meridional circulation. We use a for- 
mulation based on energy conservation along isobars similar 
to Macdcr fc Zahn (1998). In spherical polar coordinates the 
circulation takes the form 



U = U{r)P2{cos9)er + V{r)^^^^^eo, 

cLu 

where U and V are linked by continuity so that 

V=-^Upr'U) 
6pr dr 



(9) 



(10) 
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and P2{x) 
and 



^(Sx^ — 1) is the second Legendre polynomial, 



U = Co 



m^ftgcB CppT Vad - V + V 



e,n 2-kGp J 



3Gm 



where rrica = m ll 



(11) 



the total 



27rGp J' ^ ^ "^""'^ ' -^grav, 

local = L/m, Cp is the specific heat 

capacity at constant pressure, V is the radiative temperature 
gradient, Vad is the adiabatic temperature gradient, is 
the mean molecular weight gradient and Co is a calibration 
constant. 

Differential rotation arises in stars because of hydro- 
static structural evolution, mass loss and meridional cir- 
culation. This leads to hydrodynamic instabilities that re- 
distribute angular momentum. The resulting turbulence is 
much stronger horizontally than vertically and so variables 
are assumed to be roughly constant over isobars. Specifi- 
cally we can describe the angular velocity distribution by 
n = n{r). 

Taking into account all of the processes described in 
section [27T] we use, a s the evolution equation for the angular 
velocitv (|Za"hnlll992h , 



dt 



5r'' 



d{r^nu) 1_ 



dr 



d_ 

dr 



dr 



D. 



shear ^ 



and for the chemical evolution 



dcj 
dt 



J.2 



(-Dshcar + -Dcfi + Dn=o) r 



4.8^ 

dr 
idn 
dr 



dr 



(12) 



(13) 



where Ci is the abundance of element i. The diffusion coeffi- 
cient Dconv is non-zero only in convective zones and -Dshoar 
and Deff are non-zero only in radiative zones. The coefficient 
feff describes the effective diffusion of chemical elements 
because of the interaction between horizontal diffusion and 
meridional circulation. The variables U , Dshoar and D^s dif- 
fer between the two test cases as described in section [2^ 



2.2 Star maker 

STARMAKER is a p opulation synthesis code described by 
iBrott et al.l(|2011lj ). It was originally des i gned to work with 
the evolutionary models of iBrott et al.l (|2011ar ) . We have 
adapted it for use with ROSE stellar evolution models. Based 
on a grid of evolutionary models, STARMAKER interpolates 
for stellar properties given an initial mass, initial surface 
velocity and age. These are chosen at random according to 
user-defined distribution functions. Each simulated star is 
assigned a random orientation in space. The newly gener- 
ated sample can subsequently be filtered according to ob- 
servational selection effects to enable comparison with ob- 
served samples. Such effects are not applied in this paper 
because the differences we describe are not strongly affected 
when they are. 

In this study we limit ourselves to models on the main 



sequence, so if an age beyond the main-sequence is assigned 
to the model it is excluded from the simulation. Our initial 
masses are distributed with a Salpeter initial mass function. 
The initial su rface rotation velocity distribution is that of 
iDufton et al. I I2OO6.) for Galactic B-type stars, a Gaussian 
function truncated at zero with mean p, = 175km s~^ and 
standard deviation a = 94km s~^. 



2.3 Test cases 

We consider two models for comparison (cases 1 and 2 of 
Paper I). For both models we evolve a grid of stars with 
masses between 3 and 100 Mq and initial equatorial surface 
rotation velocities between and 600kms~^. The zero age 
main sequence is the point of minimum luminosity at the 
onset of hydrogen burning. The masses computed are 



m/Me € {3, 4, 5, 6, 7, 8, 9, 10, 12, 15, 20, 25, 

30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 

85, 90, 95, 100} (14) 

and for each mass the initial surface velocities used are 



t)ini/km s~^ G {0, 50, 100, 150, 200, 250, 300, 
350,400,450,500,550, 600}, 



(15) 



except when the rotation velocity would be too close to criti- 
cal rotation to achieve numerical convergence. This becomes 
more difficult for stars less massive than 10 Mq . Convergence 
can be achieved for a 10 Mq star rotating faster than 95% of 
critical rotation, although the assumptions of the model are 
likely to become invalid this close to critical. For a 3 Mq star 
the limit for convergence is close to 70% of critical rotation. 
Both the case 1 and case 2 models for each mass and initial 
surface velocity must reach the end of the main sequence 
for either of them to be used in the grid. The end of the 
main sequence is the point of maximum temperature before 
a star moves onto the Hertzsprung gap. Each model evolved 
is plotted in Fig. [T] For both models, the diffusion of angular 
momentum in convective zones is determined by the char- 
acteristic eddy viscosity given by mixing length theory such 
that Dconv = -Dniit = |t^mit^mit. The position of the convec- 
tive boundary is determined by the Schwarzschild criterion 
and although the code includes a model for convective over- 
shooting, we do not use it in this paper. Unlike for Paper 1, 
we only consider the case in which the convective core tends 
to a state of solid body rotation. In section 13.51 we exam- 
ine the effects of changing the free parameters associated 
with the model. Models generated with this calibration are 
referred to as case 2b. In this paper we generate models 
with two different metallicities. Galactic and Large Magel- 
lanic Cloud (LMC), as defined bv lBrott et"al] (|2011al '). Other 
than the initial composition, the input physics is the same 
for both metallicities. However, for clarity, we distinguish 
models that use LMC metallicity by referring to them with 
a superscript 'Z' (e.g. case 1^). 
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Agc/yr 


Number of 
excluded stars 


Maximum mass 
(case 1)/M© 


Maximum mass 
(case 2)/M0 


5 X 10''' 


250625 


42.2 


41.6 


10^ 


831418 


21.6 


22.3 


2 X 10^ 


1681936 


14.1 


13.7 


5 X 10^ 


3550414 


8.2 


7.4 



Table 1. The properties of different single-aged stellar popula- 
tions used in section[3] The original size of the population in each 
case is lO'^ stars. Each population is generated by an instanta- 
neous burst of star formation at t = 0. The first column shows 
the age of the simulated population. The second column shows 
the number of stars that have reached the end of the main se- 
quence and so are excluded from the sample. The third and fourth 
columns show the mass of the most massive star remaining in the 
sample at the given age for case 1 and case 2 respectively. 



Figure 1. Grid of initial models, in initial mass— initial equatorial 
velocity space, used for simulating stellar populations. 



and 



2.3.1 case 1 

Our case 1 mod el uses the formulation for Dshoar of 
iTalon et all l|l997h . 



^shcar — Cq 



where 



and 



Hp 



N?, = ^'^^ 



dlnp 
dlnT 



(Vad - V) 



Hp 



d\np\ dlnfi 



(16) 



(17) 



(18) 



^ 9 In /i Jpj, dlnP' 

As for paper 1 we follow lMaede3 (|2003h by taking the critical 
Richardson number, Ric = (0.8836)^/2. We have also cho- 
sen Co so that we reproduce the terminal-age main-sequence 
(TAMS) nitrogen enrichment of a 40 Mq star initially ro- 
tating at 2 7 0kms~ ^ with Galactic composition given by 
iBrott et al.l |201l3)- The effective diffusion coefficient Dcfi 



Deft 



and we take 



where 



307?h' 

Dh = 0.134r {rnV [2V ~ aU])^ 



(19) 
(20) 
(21) 



2.3.2 case 2 

Our case 2 model is that of iHeger et al] l|2000l '). In this case 
[7 = because circulation is treated as a purely diffusive 
process. The details of the various diffusion coefficients are 
extensive so we refer the reader to the original paper. With 
their notation the diffusion coefficients are 



Acm + -Ddsi + Dsni + -Dssi + Des + Dgsf (22) 



D,a = (/c - 1)(Ddsi + Dsni + Dssi + Des + I>gsf), (23) 

where ea ch Dj corresponds t o a different hydrodynamical in- 
stability. HegereFal] ([ioO^) take /c = 1/30 and we use this 
too. We also use ffi = 0. The consequences of this are dis- 
cussed in Paper 1. Unless otherwise stated, we calibrate this 
model by scaling Des, the dominant diffusion coefficient, so 
that the nitrogen enrichment of a 20 M©, solar metallicity 
star with initial surface angular velocity of 300 kms~^ is the 
same as for case 1 at the terminal-age main sequence. 



2.4 Stellar populations 

Throughout this paper we use a variety of populations at 
different ages with different star formation histories. The 
main reason for this is that, as a population ages, the mass 
of the most massive stars remaining in the main-sequence 
population decreases, whilst stars much less massive than 
the maximum mass have not had sufficient time to produce 
significant nitrogen enrichment. The combination of these 
tendencies allows us to follow how the amount of enrich- 
ment varies with mass. This applies specifically to clusters 
in which we expect the range of ages of the stars to be small 
compared to the age of the cluster. It is important to note 
that rotation can significantly affect the upper bound to the 
mass of stars in the population. The populations used in this 
paper are listed in table ([T]). For the figures in section [3] we 
compare the data using 2D histograms for which we have 
separated the data into a grid of 50 x 50 bins. The number 
of stars in each bin divided by the total number of stars in 
case 1 is ni and similarly 712 for case 2. In each comparison 
we reduce the size of the larger population to be the same 
size as the other by randomly removing stars. 



3 RESULTS 

We simulated stellar populations at a number of ages and 
metallicities for each model and found a number of signifi- 
cant differences. 
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Figure 2. Hcrtzsprung-Russell diagrams for a population of stars at age 2 X 10^ yr. In each case we have plotted the evolution of 5 Mq 
and 10 Mq stars initially rotating at 300 km s^^. There is some slight variation between the two cases at the end of the main sequence 
(highest luminosity) but the effect is small. Otherwise there is no obvious difference between the results produced in cases 1 and 2. When 
comparing the two populations, the addition of lO"'^ in the denominator is to avoid division by zero in unpopulated bins. 



3.1 The HR diagram 

When we look at the effect of the two models of rotation 
on stars in the HR diagram we find that there is very lit- 
tle difference between them. Whilst there is variation in the 
TAMS temperature and luminosities of the stars in each 
case, the difference is small and, for a single burst of star 
formation, only affects a handful of stars in the population 
at any given time. For most of a star's lifetime the predicted 
position in the HR diagram is sufficiently similar between 
the two cases that the difference in the population cannot 
be distinguished. Fig. [2] shows the HR diagrams for cases 1 
and 2 for simulated clusters with an age of 2 x lO'^ yr. Apart 
from slightly different degrees of broadening at the main- 
sequence turn off, there is no difference between the two 
cases. This is true at all ages and if we simulate a popu- 
lation of stars with continuous star formation we still find 
only slight distinctions between the two cases. This doesn't 
mean that mass determinations of rotating stars from their 
surface rotation, temperatures and luminosities are unaf- 
fected by the specific physics of the model. For individual 
models, the difference can be significant but the cumulative 
effect has little impact on the population as a whole. It is 
also important to note that the maximum mass of stars re- 
maining in the sample varies between the cases because the 
main-sequence lifetimes are different. This means that care 
must be taken when identifying a cluster's age with respect 
to its most massive members if the cluster contains rapid 
rotators. 



3.2 Velocity distribution evolution 

Because of variations in the amount of mixing and the evo- 
lutionary timescale between the two cases, we might expect 
differences between the distribution of rotation rates as the 
populations evolve. In Fig. [3] we plot the velocity distribu- 
tion of the remaining stars in the single-aged populations at 
2 X lO^yr. This is the typical shape of the distribution at all 
ages considered and we see that there is very little difference 
between the two cases. 



3.3 The Hunter diagram 

The effect of rotation on the surface abundance of various 
isotopes is an extremely important tracer of the effects of ro- 
tation. In section lT^ we showed that, using alternate models 
of rotation, we find only small effects on the velocity distri- 
bution in stellar populations. We now consider the effect 
of rotation on the surface abundance of nitrogen. We could 
make similar conclusions about other chemical elements but 
their usefulness depends on the accuracy to which they can 
be measured and the availability of data. For example in Pa- 
per 1 we discussed the effect of rotation on the surface abun- 
dance of helium-3 but this is difficult to measure and so is 
not pa rticu larly useful i n this d iscussion. iFrischknecht et al] 
(|2010l ) and iBrott et al.l l|2011al ) also consider how rotation 
is likely to affect the surface abundances of light elements. 
If we look at a plot of the surface abundance against sur- 
face r otation rate, commo nly referred to as the Hunter dia- 
gram l|Hunter et al.ll2009l ). for different ages (Fig. 2} we see 
that there are some very clear differences between the two 
cases. At each age, the most massive stars remaining in the 
population dominate the enriched stars. Stars more massive 
than this have already evolved off the main sequence. The 
less massive stars in the population evolve more slowly and 
so have not had enough time to become enriched. At early 
times the populations are very similar except that case 1 
predicts rather more enrichment for stars rotating slower 
than 200kms~^ while case 2 predicts more enrichment of 
the most rapidly rotating stars. As the population ages, the 
amount of enrichment in case 1 stays roughly the same but 
the amount of enrichment in case 2 drops off slowly followed 
by a large drop between 2 x lO'^yr and 5 x lO'^yr. It is here, 
where only stars less massive than 8.2 Mq remain, that the 
difference between the models is clearest. However, even at 
2 X lO'^yr, we can see that there are far more enriched stars 
in case 1 than case 2 compared with earlier times. 

If we consider the case dependence of the Hunter dia- 
grams for a population of stars with a continuous star for- 
mation history we find much the same thing. Fig. \5\ shows 
that, although both cases follow a trend of more frequent en- 
richment in stars with higher rotation rates, case 1 produces 
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Figure 3. Velocity distribution of stars in cases 1 and 2 both at an age of 2 X lO^yr for the populations described in table [T] There is 
no marked difference between the two distributions. 



far more enriched stars than case 2. This is not surprising 
because we found far less mixing in case-2 stars at lower 
masses than in case 1 and it is these stars that dominate the 
population. We could have instead chosen to calibrate case 2 
so that there were more mixing in low-mass stars but this 
would inevitably lead to a worse match in the populations 
elsewhere, perhaps in the enrichment of rapidly rotating very 
massive stars (M > 40 M0). We discuss this further in sec- 
tion 13.51 Also important to consider is the effect of metal- 
licity. In Paper 1 we found a reversal in the trend of less 
mixing in case-2, low- mass stars (M < 20 Mq). Increasing 
the mixing here to bring the two populations in line would 
make the low-metallicity agreement far worse. We discuss 
this in section [T6l 



3.4 Effective surface gravity and enrichment 

As we suggested in Paper 1, the difference between the 
two cases can be seen most clearly by considering differ- 
ent masses of stars. In our discussion of the Hunter diagram 
we have appealed to the single-aged population of stars to 
differentiate between stellar masses as the population ages. 
Unfortunately, determining the mass of rotating stars self- 
consistently is difficult because of the degeneracy that arises 
owing to rotation. Fig. [6] shows the typical relationship be- 
tween mass and effective surface gravity in a simulated pop- 
ulation. There is a strong correlation between the two but 
rotation causes degeneracy so estimates of the mass from 
effective gravity alone could be wrong by up to 7 Mq in this 
case. The correlation does not persist in the case of continu- 



ous star formation. Use of the effective surface gravity is also 
advantageous because it can be directly determined spectro- 
scopically. However, caution is necessary for rapid rotators 
because the effective gravity is not uniform across the stel- 
lar surface ijVon Zeipe]||l924l ). The Hunter diagram suffers 
from the problem that, even for simple stellar populations 
like this one, stars exist in all regions of the diagram and 
the population has few clear boundaries. If we look at the 
variation of effective surface gravity with nitrogen enrich- 
ment the difference between the models becomes very clear 
(Fig. [7]). There are sharp curves that bound the upper and 
lower effective surface gravities of the population. The lower 
bound occurs because stars evolve rapidly into giants with 
much lower surface gravity after this limit. The upper bound 
occurs because younger stars with higher surface gravities 
haven't evolved to the point where their surface nitrogen is 
enriched. There are features which distinguish the two popu- 
lations at each age. For young populations (5 x lO^'yr) case 2 
has a higher upper bound for nitrogen enrichment and there 
is a much broader range of surface gravities than in case 1. 
For older populations (lO^yr and 2 x lO^yr) case 2 predicts 
generally lower values for the surface gravity. Finally for old 
populations (5 x lO^yr) the difference becomes very stark. 
The amount of mixing in case 2 drops off dramatically com- 
pared to case 1 while we still predict much lower values for 
the surface gravity in rapid rotators. 

When we consider a population of stars with continuous 
star formation history, the difference in the populations is 
still clear. Interestingly, unlike the Hunter diagram, this vi- 
sualisation actually highlights the similarities as well as the 
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Figure 4. Hunter diagrams for single-aged populations of stars. Prom top to bottom, the four rows of figures correspond to 5 X lO^yr, 
10''^yr, 2 X lO'^yr and 5 X lO^yr. At early times the two cases are similar with slightly more enrichment of the fastest rotators in case 2 
and more enrichment of stars rotating more slowly than 200kms~^ in case 1. By 2 X lO^yr we see many more enriched stars in case 1 
and by 5 X lO^yr the amount of mixing in case 2 has dropped off dramatically. The jagged right hand edge of the populations is a result 
of the grid geometry and the mass-independence of the initial rotation velocity distribution. Neither affects the large difference we see in 
the populations at late times. 
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Figure 5. Hunter diagrams for a population of stars witli continuous star formation. In cacli case we have plotted the evolution of 
IOMq stars with initial surface rotation velocities 200kms~^, 300kms~^ and 400kms~^. Despite showing good agreement at low 
surface rotation rates, case 1 has many more fast-rotating highly enriched stars. However, this difference can often be accounted for by 
recalibration of the mixing coefficients and is difficult to observe owing to the rarity of rapidly rotating high-mass stars which occupy 
this region of the plot. 
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Figure 6. Simulated single-aged stellar population at lO'^yr in 
case 1. The plot shows the strong correlation between mass and 
surface gravity. This relation only holds when the population has 
a single age. Because of rotation, the relation is degenerate and a 
measurement of the surface gravity corresponds to stellar masses 
with a range of up to 7M0. 



differences between the two cases. Fig.[S]shows that stars in 
both models are confined to a similar band of effective grav- 
ities and their range of surface abundances are very similar. 
The main difference between the two cases, apart from the 
increased frequency of enriched stars in case 1 which we 
saw in section [3]3j is the confinement of the enriched case-1 
stars to a distinct band. This contrasts to case 2 for which 
the stars are spread much more evenly across their range of 
enrichment. 



3.5 Recalibration 

We have thus far described the differences that arise between 
the two test cases under a specific calibration of the mixing. 
However, within each case there is the fiexibility to calibrate 
to some degree the amount of mixing that arises because of 
rotation. We chose in our initial calibration to match the 
TAMS nitrogen enrichment of 20 Mq stars initially rotating 
at i; = 300kms~^. This is a reasonably good fit for stars of 
M > 15 Mq but for smaller masses the amount of mixing 
in case 2 drops off rapidly. Now suppose instead we had 
chosen to match the TAMS nitrogen enrichment of a star 
with initial mass M = 10 M© and v — 200kms~^. We refer 
to this model as case 2b. This is more representativ e of the 
stars observed in the VLT-FLAMES survey (Dufton et all 
l2006l ) and so should produce mixing in line with the bulk of 
the population. We discuss the VLT-FLAMES survey data 
in relation to our simulated populations in sect ion [3?7l Fig. [9] 
shows the Hunter diagram for the new sample. We see that 
the agreement is better in the Hunter diagram but case 2b 
still can't produce the tightly confined bulk of enriched stars 
seen in case 1. Also, the maximum enrichment observed in 
case 2b is now far greater than in case 1. 

Although we can't directly measure the mass of stars, 
it is instructive to examine where the main differences in 
our sample arise. Fig. [10] shows the distribution of nitrogen 
enrichment by mass in case 1, case 2 and case 2b. In the 
first instance, the agreement between the two models is rea- 
sonable for stars more massive than around 15 M© but the 
mixing in case 2 drops off rapidly for lower masses as we 
observed in section \3A\ For case 2b the agreement holds to 
much lower masses except now there are far more highly en- 
riched stars with mass 5 Mq < M < 20 Mq when compared 
to case 1. 



3.6 Effects of metallicity 

In Paper 1 we concluded that there was significant variation 
in the differences between the two cases at different metal- 
licities. We simulated a grid of models at LMC metallic- 
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Figure 7. Distribution of the surface nitrogen enrichment against effective surface gravity in single-aged stellar populations. From top 
to bottom, the four rows of figures correspond to 5 X lO^yr, lO^yr, 2 X lO^yr and 5 X lO'^yr. At early times case 2 gives a larger spread 
of effective surface gravities and higher enrichment of the fastest rotators. At later times the maximum enrichment is similar but case 2 
predicts overall lower surface gravities than case 1. Finally at late times when only stars with mass smaller than 8.2 Mq remain, case 2 
predicts far less mixing than case 1 as well as much lower surface gravity for its fastest rotators. 
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Figure 8. Surface nitrogen enrichment against effective surface gravity of simulated populations of stars with continuous star formation. 
The distributions are both confined to a narrow band and have similar ranges for enrichment, though slightly higher in case 2. However, 
case 1 produces many more enriched stars than case 2 when they are for the large part confined to a narrow band. There are some 
edge-of-grid effects that arise because the initial rotation function is mass independent and so produces many more low-mass stars close 
to their critical limit. 




-1.5 






9 


-2 








-2.5 






8.8 


-3 






8.6 


-3.5 








-4 






8.4 


-4.5 


i 


i 
+ 


8.2 


-5.5 






8 


-e 

-6.5 






7.8 


-7 










-1.5 




9 


-2 






-2.5 




8.8 


-3 




8.6 


-3.5 






-4 £ 




8.4 


-4.5 '§ 








+ 


8.2 








-5.5 




8 


-6.5 




7.8 









200 300 

L'sini/kiiiH ^ 



Figure 9. Hunter diagrams for a population of stars undergoing continuous star formation with case 2 calibrated to give the same 
TAMS nitrogen enhancement as case 1 for a star of mass 10 Mq initially spinning with v = 200kms~^ (case 2^). In each case we have 
plotted the evolution of IOMq stars with initial surface rotation velocities 200kms~^, 300kms~^ and 400kms~^. There are now more 
moderately enriched stars in case 2i, but still far fewer than in case 1 and the upper limit for enrichment in case 2^ is now far greater 
than case 1. 



ity, composition an d initial velocity distribution, as given by 
iBrott et al.l (|2011al '). This is not as low as the low-metallicity 
case described in Paper 1 but it does allow us to compare 
our results with the data from the LMC observations in 
the VLT-FLAMES survey of massive stars. We simulated a 
population in case 1^ and case 2 at this composition with 
continuous star formation history. The Hunter diagram for 
this population is shown in Fig. 1111 

We see that the qualitative distribution of stars in the 
simulated population is similar to the solar metallicity pop- 
ulations. Case 1^ produces a much more well-defined band 
of enriched stars whereas case 2^ produces fewer enriched 
stars that have a much greater spread in abundance. If we 
consider the mass-dependence of the rotational mixing we 
find a similar decline in the amount of mixing in case 2^ 
compared to case 1^ for stars less massive than 20 M0. As 
in Paper 1, we find that the amount of mixing in stars above 
this mass is higher in case 2^ than in case 1^. In fact, the 
mixing in case 1^ decreases slightly for higher-mass stars. 



This means that, as in section [3751 an increase in the mixing 
in case 2^ is unlikely to produce a better correlation between 
the two cases. However, owing to the IMF, there are many 
more stars less massive than 20 Mq in the population, case 1 
produces many more enriched stars than case 2^. 



In Fig. [TT] we have also plotted those LMC stars for 
which the ni t rogen abundances have been determ i ned b y 
iHunter et all (|2009l '). As remarked by iHunter et al.l l|2009f) . 
there are many highly enriched, slowly rotating stars that 
are not explained by either model of rotational mixing. How- 
ever, for the remainder of observed stars we see a trend of 
increasing enrichment for higher rotation rates. On initial 
inspection, case 1^ fits the VLT-FLAMES data much more 
closely than case 2^. However, selection effects are important 
and we consider these in section [3771 
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Figure 10. Distribution of stars in case 1 and case 2 for a population witli continuous star formation history. The top row shows the 
models with the calibration of case 2 to give the same TAMS nitrogen enhancement as case 1 for a star of 20 Mq initially spinning with 
V = 300kms~^. The second row shows the same but with case 2 calibrated to give the same TAMS nitrogen enrichment as as star with 
initial mass M = 10 Mq and surface rotation v = 200kms~^ (case 21,). The two cases agree for masses greater than 15 Mq for the 
original calibration. The agreement continues to lower masses for the second calibration but now there are more highly enriched stars in 
case 2b in the mass range 5Mq< M < 20 Mq. 



3.7 Selection effects in the VLT-FLAMES survey 

The most common data set used to test rotational mixing in 
massive s tars is th e VLT - FLAMES survey o f massive stars 
ijEvans et aL. .2005i . l2006l : iDufton et aDl2006l ) owing to the 
number of stars sampled and the detailed determination of 
surface composition. We repeated our population synthe- 
sis as in section 13.61 for a continuous population of LMC- 
metallicity stars but we have included the selection effects 
which affect the stars in the VLT-FLAMES survey so that 
we may compare the distribu tions more directly with those 
found byHuntereTal] l|2009l ). The selection criteria we used 
are that of th e cluster Ni l in th e LMC. For a detailed 
description see iBrott et al.l (|2011bl ). Stars are excluded if 
their visual magnitude is greater than 15.34, if they are 
hotter than 35,000 K, if their surface gravity is less than 
lO^'^cms"^ or they are rotating faster than 90% of their 
critical rotation rate. In addition, a random error is applied 
to logif)[N / H] selected randomly from a Gaussian with stan- 
dard deviation a — 0.2. 

The simulated population produced after we apply the 
selection effe cts is shown in Fig. [12] along with the LMC 
stars data of llfunter et al.l ((2009). We show the population 
in cases 1^ , 2^ and 2^ (the LMC analogue of case 2b with the 
calibration described in section [33|) . We see that, contrary to 
our discussion in sect ion [TBI the differences between the var- 



ious models are now far less apparent. Compared with the 
other two cases, case 2^ predicts many more less-enriched 
fast rotators than are observed and the amount of mixing 
is insufficient to match the observed band of enriched stars. 
Compared with section [T6] cases 1^ and 2f now both show 
a similarly good fit to the data. In both cases the band of 
predicted enriched stars is matched well by the observations. 
However, we note that we would expect to see a number of 
stars with 12 + logj^Q [A^/ff] > 8. The number of predicted 
stars in this range is greater for case 2b than case 1^. If we 
were to reduce the amount of mixing we would get too little 
nitrogen enrichment in stars with v < 200 kms^'^. This effect 
may be a result of the difficulty of measuring nitrogen abun- 
dances in this region but is otherwise difficult to resolve. A 
further increase in the mixing would exacerbate the problem 
that the upper bound to enrichment at rapid rotation is too 
high. A decrease in the mixing would mean that the band 
of enriched stars in the simulated population is less likely to 
produce a good fit to the model and would leave observed 
slowly-rotating, moderately-enriched stars that cannot be 
explained through the theoretical models. 

Unfortunately, the mass dependence of the models is 
not well refiected in the VLT-FLAMES populations. With 
the selection effects described, the masses of the LMC- 
metallicity sample are confined between 10 Mq and 20 Mq 
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Figure 11. Hunter diagrams for populations of stars undergoing continuous star formation at LMC metallicity. The two population s 
are qualitatively similar to the populations at solar metallicity. We have also plotted the LMC stars observed bv lHunter et al] l l2009t) . 
Case 1^ predicts a very confined distribution of enriched stars, whereas case 2^ predicts a much wider spread of enrichment. The stars 
at the left-hand edge of the diagram can't be explained by rotational mixing alone. 
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Figure 12. Hunter diagrams for populations of stars undergoing continuous star formation at LMC metallicity for a number of different 
models with selection effects applied as described in section PS.?! Populations have been simulated for cases fZ, 2^ an d 2g. The popu lat ions 
are qualitatively similar to the populations at solar metallicity. We have also plotted the LMC stars observed bv lHunter et al] l l2009fl . 
The slowly-rotating, highly-enriched stars at the left-hand edge of the diagram cannot be explained by rotational mixing alone. 



as shown in Fig. 1131 We find a similarly narrow range when 
we consider Galactic stars under similar selection effects. 
Therefore, any simulated population where we include the 
selection effects of VLT-FLAMES captures the rotational- 
dependence of the model but for only a small fraction of the 
mass-dependence which is where we have found the biggest 
differences between our two cases. 



4 CONCLUSIONS 

Rotation has many effects on stellar evolution. Some of 
these, such as that on the surface temperature, are be- 
cause of rotation but only vary significantly between differ- 
ent models towards the end of the main sequence. Others, 
such as that on the surface rotation velocity, may evolve dif- 
ferently throughout the main sequence according to different 
models but do not produce significant changes in the distri- 
bution of stars in a simulated population. Therefore these 
properties alone are largely unhelpful to distinguish between 
the different implementations of rotation in stellar models. 



It has been observed that the surface abundances of 
several chemical elements change significantly because of 
rotation. The degree to which this happens in the theoreti- 
cal models strongly depends on which particular model for 
stellar rotation is used and which constraints are used to 
calibrate them. Recently, the Hunter diagram has been the 
favoured diagnostic tool for analysing stellar rotation be- 
cause it shows a clear connection between the surface rota- 
tion of a s tar and its surface enrichment. Our model based 
on that of iTalon et al] (|l997l '). case 1, shows a similar order 
of magnitude enrichment for all masses. On t he other hand 
our model based on that o f lHeger et al] (|2000l ). case 2, shows 
a steep decline in the amount of enrichment around ISM©. 
We can account for this to a degree by adjusting the cali- 
bration of the models but we see that increasing the mixing 
in case 2, so that low-mass stars show similar enrichment to 
those of case 1, then leads to a much higher maximum en- 
richment for a case-2 population than a case-1 population. 
This suggests that studies should focus on stars either side 
of this mass limit. 
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Figure 13. Distribution of masses in a simulated population of 
case-1^ stars at LMC metallicity with the inclusion of selection 
effects. The distribution is very similar in each case. We see that 
the distribution peaks strongly around 12 Mq. 



Because the two models are very different for different 
mass ranges, the effective gravity is a sensible tool to investi- 
gate the mass-dependence of the mixing strength, ft is very 
difficult to self-consistently infer the mass of a star from a 
luminosity-temperature-rotation relation because any such 
relation depends on the model used for rotation. The effec- 
tive gravity however can be derived directly from spectra 
and, despite some degeneracy, there is a strong relation be- 
tween it and the stellar mass. This can be usefully applied 
to the study of rotational mixing. We have shown that there 
are very clear differences between the synthetic populations 
produced by each model when nitrogen enrichment is plot- 
ted against surface gravity. In the case of a population in 
which all the stars have the same age, they are confined to a 
very specific region. Stars with surface gravity below a cer- 
tain limit evolve into giants and beyond on a short timescale 
compared to their main-sequence lifetime. Stars with sur- 
face gravity above another limit are less massive and have 
not had long enough to become enriched. This effect persists 
even in the case of a population of stars with continuous star 
formation. Case 1 predicts that nitrogen is enriched over a 
narrower range of surface gravities than case 2. We see that 
case 1 produces many more moderately enriched stars than 
case 2 but the maximum enrichment in case 1 is lower than 
in case 2. 

We have also shown that similar trends appear at dif- 
ferent metallicities. We simulated populations in both cases 
with continuous star formation history at LMC metallicity. 
The qualitative distribution of stars in each case was sim- 
ilar to that at Galactic metallicity. Case 1^ still produced 
a confined band of enriched stars in the Hunter diagram, 
whereas case 2^ produced a much greater spread. Case 1^ 
also produced many more highly-enriched stars than case 2^ 
although the maximum enrichment in case 2^ was higher 
than in case 1^. Similarly to solar metallicity, this is be- 
cause of the mass dependence of the mixing in each case. 



The decline in the amount of mixing in case 2 begins at 
higher masses at low metallicity (20 Mq at LMC rather than 
15 Mq at solar metallicity) and the mixing in stars above 
this mass is relatively constant for stars in case 2^, whereas 
case-1^ stars show slightly less enrichment as the mass of 
stars increases. This may not be indicative of the strength 
of the rotational mixing but rather that, because the main- 
sequence lifetime of the stars decreases with increasing mass, 
there is less time to transport nitrogen from the core to the 
surface. 



When we compare the simulated populations to the 
LMC-metallicity stars in the VLT-FLAMES survey we find 
that both cases 1^ and 2^, which uses the second calibra- 
tion for case 2 described in section 13.51 give a reasonable fit 
to the observed data and it is difficult to determine which 
fits the data more closely. Unfortunately the range of initial 
masses that remain in the simulated samples after we apply 
the selection effects is extremely narrow, between 10 Mq and 
20 Mq in the LMC populations. This means that it is diffi- 
cult to observe the large difference in the mass-dependence 
of the models. Hence, a close fit between the Hunter dia- 
grams for a simulated population and the data of the VLT- 
FLAMES survey is a useful test of models for rotational 
mixing but cannot establish the validity of a model by itself. 
As we have seen, models for rotational mixing with very dif- 
ferent mass-dependencies can reproduce similarly good fits 
to the current data. Whilst, in the future, determinations 
of the nitrogen abundance in a wider mass-range of stars 
may help solve this problem, it is likely that examining the 
enrichment and d epletion of o t her ele ments will be neces- 
sary. For example, iBrott et al.l l|2011al ) look at the effect of 
rotation on the surface Boron abundance. Different initial 
rotation velocity distributions may also have a significant 
effect on the simulated populations. For example an initial 
distribution which is a function of Q/ilciit would prevent the 
over-abundance of rapid rotators at low masses. 



Despite sharing many similar features, it is unreason- 
able to expect that two different models for stellar rotation 
can produce identical qualitative results for an extended 
range of masses, rotation rates and metallicities. We have 
shown that, whilst the two models agree for stars more 
massive than 15 Mq , there is much less agreement for less 
massive stars. Furthermore, we have only thus far made 
a comparison of two particular models for stellar rotation 
chosen from the many available. In particular we haven't 
yet in cluded mode ls based around the Taylo r-Spruit dy- 
namo jSpruitll2002l ) such as that investigated b vlBrott et al 
dSOljj). Whilst this is a similar model to that of lHeger et al 
(2000) it produces very different results owing to the in- 
clusion of magnetic fields. This extension is supported by 



our consideration of different metallicity regimes. Whilst the 
better fit of case 1 suggests that models in which meridional 
circulation is treated advectively and diffusion comes solely 
from h ydrodynamical insta bilities are more realistic, the re- 
sults of iBrott et al.l l|2011b[ ) also produce a reasonable fit to 
observed LMC stars. To distinguish between these models, 
and others, it is necessary to continue with this analysis and 
extend it to different masses and metallicities as more data 
becomes available. 
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